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ABSTRACT 

The acceleration time scale for the process of first-order Fermi acceleration in rela- 
tivistic shock waves with oblique magnetic field configurations is investigated by the 
method of Monte Carlo particle simulations. We discuss the differences in derivation 
of the cosmic ray acceleration time scale for non-relativistic and relativistic shocks. We 
demonstrate the presence of correlation between the particle energy gain at interaction 
with the shock and the respective time elapsed since the previous interaction. Because 
of that any derivation of the acceleration time scale can not use the distribution of 
energy gains and the distribution of times separately. The time scale discussed in the 

(c) 

present paper, Tqcc, is the one describing the rate of change of the particle spectrum 
cut-off energy in the time dependent evolution. It is derived using a simplified method 
involving small amplitude particle momentum scattering and intended to model the 
situations with anisotropic cosmic ray distributions. We consider shocks with parallel, 
as well as oblique, sub- and super-luminal magnetic field configurations with finite 
amplitude perturbations, SB. At parallel shocks Tacl diminishes with the growing 
perturbation amplitude and shock velocity Ui . Another feature discovered in oblique 

shocks are non-monotonic changes oiTacl with 6B. The effect arises due to the particle 
cross-field diffusion. The acceleration process leading to power-law spectra is possible 

fc) 

in super-luminal shocks only in the presence of large amplitude turbulence. Then, Tacc 
always increases with increasing SB. In some of the considered shocks the accelera- 
tion time scale can be shorter than the particle gyroperiod upstream the shock. We 
also indicate the relation existing for relativistic shocks between the acceleration time 
scale and the particle spectral index. A short discussion of the numerical approach 
modelling the pitch angle diffusion versus the large angle momentum scattering is 
given. We stress the importance of the proper evaluation of the effective magnetic 
field (including the perturbed component) in simulations involving discrete particle 
momentum scattering. 

Key words: cosmic rays - shock waves - acceleration time scale - Fermi acceleration 



by Kirk & Schneider (1987a; see also Kirk 1988). By extend- 
ing the diffusion approximation to higher order terms in the 
anisotropy of particle distribution, they obtained solutions 
to a kinetic equation of the Fokker-Planck type with the 
isotropic form of pitch angle diffusion coefficient. Next, Kirk 
& Schneider (1988) extended the analysis by involving both 
diffusion and large-angle scattering in particle pitch angle. 
They discovered that - in relativistic shock waves - the pres- 
ence of scattering can substantially modify the spectrum of 
accelerated particles. An extension of Kirk & Schneider's 



A consistent method to tackle the problem of first-order 
Fermi acceleration in relativistic shock waves was conceived 



1 INTRODUCTION 



(1987a) approach to more general conditions in the shock 
was given by Heavens & Drury (1988), who took into consid- 
eration the fluid dynamics of relativistic shock waves. They 
also noted that the resulting particle spectral indices depend 
on the perturbations spectrum near the shock, in contrast 
to the non-relativistic case. Kirk & Heavens (1989) consid- 
ered the acceleration process in shocks with magnetic fields 
oblique to the shock normal (see also Ballard & Heavens 
1991). They showed, contrary to the non-relativistic results 
again, that such shocks led to flatter spectra than the par- 
allel ones. Their work relied on the assumption of adiabatic 
invariant p\/B conservation for particles interacting with 
the shock, which restricted the considerations to the case of 
nearly uniform magnetic fields upstream and downstream of 
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the shock. A different approach to particle acceleration was 
presented by Begelman & Kirk (1990), who noted that in 
relativistic shocks most field configurations lead to super- 
luminal conditions for the acceleration process. In such con- 
ditions, particles are accelerated in a single shock transmis- 
sion by drifting parallel to the electric field present in the 
shock. Begelman & Kirk showed that there is more efficient 
acceleration in relativistic conditions than that predicted by 
the simple adiabatic theory. The acceleration process in the 
presence of finite amplitude perturbations of the magnetic 
field was considered by Ostrowski (1991; 1993) and Ballard 
& Heavens (1992). Ostrowski considered a particle accelera- 
tion process in the relativistic shocks with oblique magnetic 
fields in the presence of field perturbations, where the as- 
sumption p\/ B = const was no longer valid. To derive parti- 
cle spectral indices he used a method of particle Monte Carlo 
simulations and noted that the spectral index was not a 
monotonic function of the perturbation amplitude, enabling 
the occurrence of steeper spectra than those for the limits of 
small and large perturbations. It was also revealed that the 
conditions leading to very flat spectra involve an energetic 
particle density jump at the shock. The acceleration process 
in the case of a perpendicular shock shows a transition be- 
tween the compressive acceleration described by Begelman 
& Kirk (1990) and, for larger perturbations, the regime al- 
lowing for formation of the wide range power-law spectrum. 
The analogous simulations by Ballard & Heavens (1992) for 
highly disordered background magnetic fields show system- 
atically steeper spectra in comparison to the above results, 
as discussed by Ostrowski (1993). The particle spectrum for- 
mation in the presence of non-linear coupling of accelerated 
particles to the plasma flow has been commented by Os- 
trowski (1994a). The review of the above work is presented 
by Ostrowski (1994b). 

The shock waves propagating with relativistic velocities 
rise also interesting questions concerning to the cosmic ray 
acceleration time scale, Tacc- To date, however, there is only 
somewhat superficial information available on that problem. 
A simple comparison to the non-relativistic formula based on 
numerical simulations shows that Tacc relatively decreases 
with increasing shock velocity for parallel (Quenby & Lieu 
1989; EUison et al. 1990) and oblique (Takahara & Tera- 
sawa 1990; Newman et al. 1992; Drolias & Quenby 1994; 
Lieu et al. 1994; Quenby & Drolias 1995; Naito & Takahara 
1995) shocks. However, the numerical approaches used there, 
based on assuming the particle isotropization at each scat- 
tering, neglect or underestimate a rather significant factor 
controlling the acceleration process - the particle anisotropy. 
Elhson et al. (1990) and Naito & Takahara (1995) included 
also derivations applying the pitch angle diffusion approach. 
The calculations of Ellison et al. for parallel shocks show 
similar results to the ones they obtained with large ampli- 
tude scattering. In the shock with velocity 0.98 c the accel- 
eration time scale is redued on the factor ~ 3 with respect 
of the non-relativistic formula. Naito & Takahara considered 
shocks with oblique magnetic fields. They confirmed reduc- 
tion of the acceleration time scale with increasing inclina- 
tion of the magnetic field, derived earlier for non-relativistic 
shocks (Ostrowski 1988). However their approach neglected 
effects of particle cross field diffusion and assumed the adi- 
abatic invariant conservation at particle interactions with 
the shock. These two simplifications limit their results to 



the cases with small amplitude turbulence near the shocku. 
One should also note that comparing in some of the men- 
tioned papers the derived time scales to the non-relativistic 
expression does not have any clear physical meaning when 
deahng with relativistic shocks. 

In order to consider the role of particle anisotropic dis- 
tributions and different configurations of the magnetic field 
in shocks the present work is based on the small angle parti- 
cle momentum scattering approach described by Ostrowski 
(1991). It enables to model the effects of cross-field diffu- 
sion, important in shocks with oblique magnetic fields. Let 
us note (cf. Ostrowski 1993) that this code allows for a rea- 
sonable description of particle transport in the presence of 
large amplitude magnetic field perturbations also. In Sec- 
tion 2 below, we discuss the differences in derivation of the 
cosmic ray acceleration time scale for non-relativistic and 
relativistic shock waves. We demonstrate the existence of 
noticeable correlations of the particle energy gains at inter- 
actions with the shock and the respective times elapsed since 
the previous interaction. Because of that any derivation of 
the acceleration time scale cannot use the distribution of 
energy gains and the distribution of times separately. We 
define the acceleration time scale rice, as the one describing 
the rate of change of the cut-off energy in the time depen- 
dent particle spectrum evolution. Then, in Section 3, the 
performed simulations are described. We use the simplified 
method involving small amplitude particle momentum scat- 
tering and devoted to model the situations with anisotropic 
cosmic ray distributions (Kirk & Schneider 1987b, Ellison 
et al. 1990, Ostrowski 1991). In Section 4 the results are 
presented for shock waves with parallel and oblique (both, 
sub- and super-luminal) magnetic field configurations. We 
consider field perturbations with amplitudes ranging from 
very small ones up to 5B ~ B. In parallel shocks Tic' di- 
minishes with the growing perturbation amplitude and the 
shock velocity Ui. However, it is approximately constant at 
a given value oi Ui if we use the formal diffusive time scale 
as the time unit. Another qualitative feature discovered in 
oblique shocks is that due to the cross-field diffusion Tacc 
can change with 5B in a non-monotonic way. The accelera- 
tion process leading to the power-law spectrum is possible in 
super-luminal shocks only in the presence of large amplitude 
turbulence. Then, in contrast to the quasi-parallel shocks, 
rice increases with the increasing SB. In some cases with 
the oblique magnetic field configuration one may note a pos- 
sibility to have an extremely short acceleration time scales, 
comparable or even smaller than the particle gyroperiod in 
the magnetic field upstream the shock. We also demonstrate 
the existence of the coupling between the acceleration time 
scale and the particle spectral index. A form of the involved 
relation much depends on the magnetic field configuration. 
In order to evaluate some earlier simulations applying the 
large angle scattering ('LAS') model and/or the pitch angle 
diffusion ('PAD') model a short discussion of the results ob- 
tained within these two approaches is presented. We stress 



* One should note that the spatial distributions near the shock 
derived by these authors (their figures 1 and 2) do not show the 
particle density jump proved to exist in obUque relativistic shocks 
by Ostrowski (1991). It is also implcitly present in analytic deriva- 
tions of Kirk & Heavens (1989). 
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the importance of the proper evaluation of the effective mag- 
netic field (including the perturbed component) in simula- 
tions involving discrete particle momentum scattering. In 
the final Section 5 we provide a short summary and discuss 
some astrophysical consequences of our results. Some pre- 
liminary results of this work were presented in Ostrowski & 
Bednarz (1995). 

Below, the light velocity is used as the velocity unit, 
c = 1. As the considered particles are ultrarelativistic ones, 
p = E, wc often put the particle momentum for its energy. 
In the shock we label all upstream (downstream) quantities 
with the index '1' ('2'). If not otherwise indicated, the quan- 
tities arc given in their respective plasma rest frames. The 
shock normal rest frame is the one with the plasma velocity 
normal to the shock, both upstream and downstream the 
shock (cf. Begelman & Kirk 1990). In the present paper the 
acceleration time scales are always given in this particular 
frame. 



2 THE ACCELERATION TIME SCALES IN 
NON-RELATIVISTIC VERSUS 
RELATIVISTIC SHOCK WAVES 

In the case of non-relativistic shock wave, with velocity 
Ui ^ 1, the acceleration time scale can be defined as 

"""^ = ^ 

At 

where AE is the mean energy gain at particle interaction 
with the shock and At is the mean time between succes- 
sive interactions. One can use mean values hero because any 
substantial increase of particle momentum requires a large 
number of shock-particle interactions and the successive in- 
teractions are only very weakly correlated with each other. 
The respective expression for Tacc in parallel shocks. 



(2.1) 



where Ki is the respective particle spatial diffusion coeffi- 
cient, has been discussed by Lagage & Cesarsky (1983). Os- 
trowski (1988) provided the analogous scale for shocks with 
oblique magnetic fields and small amplitude magnetic field 
perturbations. It can be written in the form 



T = 



Kn,2 



U1-U2 



Ui 



U2 



Kn,2 



,2 COS-' l/'2 



(2.2) 

where the index n denotes quantities normal to the shock, 
the index || those parallel to the magnetic field, %l) is an an- 
gle between the magnetic field and the shock normal and 
Ui I cos 7/11 <C c is assumed. The terms ^ Knl{K\\ cos^ 1/)) rep- 
resent a ratio of the mean normal velocity of a particle to 
such velocity in the absence of cross-field diffusion. One may 
note that for negligible cross-field diffusion the expression 
(2.2) coincides with (2.1) if wo put K„,i for Ki (i — 1, 2). 
The case of oblique shock with finite amplitude field per- 
turbations has not been adequately discussed yet, but we 
expect the respective acceleration scale to be between the 



values given by the above formulae for Tj^cc a-nd T^cc- The 
infiuence of the particle escape boundary on the accelera- 
tion time scale and the particle spectrum is discussed by 
Ostrowski & Schlickeiser (1996). 

If the shock velocity becomes relativistic, the particle 
energy change at a single interaction with the shock can be 
comparable, or even larger than the original energy. More- 
over, after interaction with the shock, the upstream particles 
with small initial angles between its momenta and the mean 
magnetic field have a larger chance to travel far away from 
the shock. On average, such particles spend longer times 
and are able to change its pitch angles substantially until 
the next hits at the shock. Then, larger pitch angles allow 
for particle reficctions with large energy gains or transmis- 
sions downstream (cf. Ostrowski 1994b, also Ostrowski 1991, 
Lucek & Bell 1994). Therefore, correlations of the times be- 
tween successive interactions, Atdiff, the energy gains at 
these interactions, AE, and, possibly, the probability of par- 
ticle escape occur. As an example, in Fig. 1 we map the 
number of particle interaction with the shock in coordinates 
{Atdiff, AE) (a more detailed discussion of the existing cor- 
relations will be presented separately (Bednarz & Ostrowski, 
in preparation)). A cut of the presented surface at any par- 
ticular value of Atdi / / gives the distribution of energy gains 
for particles who have spent this time since the last inter- 
action with the shock. A general trend seen on the map for 
increasing Atdiff is growing the value of AE/E for the dis- 
tribution maximum. Because of that, we propose a different 
approach to derivation of the acceleration time scale with 
respect to that used for non-relativistic shocks. Usually the 
acceleration time scale is applied for derivation of the high- 
est energies occurring in the particle spectrum, characterized 
by its cut-off energy, Ec- Thus wo use this energy scale to 
define the acceleration time scale as: 



j.(c) ^ 



-Be 



(2.3) 



where Ec = dEc/dt. The rate of the cut-off energy increase 
is a well-defined quantity and the time scale (2.3) has a 
clear physical interpretation. The above definition does not 
require any limit for the energy gains of the individual parti- 
cles and all possible correlations are automatically included 
here. Prom the meaning of the definition (2.3) it follows that 
Tacl is somewhat shorter than the respective scale at the 
same energy for later times, required for the respective part 
of the spectrum to become a pure power-law (cf. Ostrowski 
& Schlickeiser 1996). One should also note that in relativis- 
tic shocks the time scale depends on the reference frame we 
use for its measurement. In the present paper the accelera- 
tion time scales are given in the respective normal shock rest 
frame. However, the applied time units re,i/c (see below) are 
defined with the use of the upstream gyration time. 



3 SIMULATIONS 

We derive estimates of Ticc basing on the Monte Carlo par- 
ticle simulations involving the small amplitude momentum 
scattering model of Ostrowski (1991) and 2500 particles in 
each run of the code. In the present considerations we dis- 
cuss the role of the mean magnetic field configuration and 
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Figure 1. Distribution of particle-shoclc interaction events for 

upstream particles, N(At^iff,AE/E), in function of the up- 
stream diffusive times, At^iff, and tlie respective energy clianges, 
AE/E. The correlation between At^aff and AE/E is repre- 
sented by a regular drift of the distribution maximum toward 
higher AE/E when increasing At^iff. An example for the paral- 
lel shock with U\ = 0.5, Vi = 1° and weak scattering conditions 
(kx/k|| = 1-6 ■ 10~®) is presented. 



the amount of particle scattering. In order to avoid the ad- 
ditional ofTccts of the varying shock compression due to the 
presence of different magnetic field configurations we take 
the field as a trace one, without any dynamical effects on the 
plasma flow. The shock compression, as seen in the normal 
shock rest frame, r, is derived from the approximate formu- 
lae presented by Heavens & Drury (1989). For illustration 
of the results, in the present paper we consider the shock 
waves propagating in the cold electron-proton plasma. For 
the mean magnetic field Bi taken in the upstream plasma 
rest frame and inclined at the angle with respect to the 
shock normal wc derive its downstream value and inclina- 
tion, B2 and 1/^2, with the use of jump conditions presented 
for relativistic shocks by e.g. Appl & Camenzind (1988): 

B2 = Bi Vcos2 Vi + sin^ V^i , (3.1) 

tan ip2 = R tan ipi , (3-2) 

where 7? = r 71/72 and the Lorentz factors 7^ = 1/ y^l-~Ui 
(i = 1, 2). These formulae are valid for both sub- and super- 
luminal magnetic field configurations. 

We model particle trajectory perturbations by intro- 
ducing small-angle random momentum scattering along the 
mean-field trajectory. The particle momentum scattering 
distribution is uniform within a cone wide at AQ (<< 1), 
along the original momentum direction. The present simu- 
lations use a constant value of An = 0.173 (= 10°). Scat- 
tering events are at discrete instants, equally spaced in time 
as measured in the units of the respective rg,i/c {i = 1, 2). 
Here we affix a gyroradius with the index 'g' when it is a 
value given for the local uniform magnetic field component. 
Index 'e' means the effective field including the field pertur- 



bations (see below). The increasing perturbation amplitude 
is reproduced in simulations by decreasing the time period 
At between the successive scatterings. For simplicity, we use 
the same scattering pattern {Ail and At in units of rg/c) 
upstream and downstream the shock, leading to the same 
values of k±/k\\ in these regions (see, however, Ostrowski 
1993). One should note that the particle momentum scatter- 
ing due to presence of turbulent magnetic field is equivalent 
to the effective magnetic field larger than the respective uni- 
form mean component, Bi or B2. In our model, the effective 
field can be estimated £is 




It is the lower limit for actual field since the amount of 
power in perturbations with wave-lengths smaller than c At 
cannot be considered within such a simple model. Below, 
the derived acceleration time scales are presented in units of 
the formal diffusive scale To = 4(Kn,i/(7i + hin,2/U2)/c or in 
units of Te.i/c , in the shock normal rest frame. 

Our numerical calculations involve particles with mo- 
menta systematically increasing over several orders of mag- 
nitude. In order to avoid any energy dependent systematic 
effect we consider the situation with all spatial and time 
scales - defined by the diffusion coefficient, the mean time 
between scatterings and the shock velocity - proportional to 
the particle gyro-radius, rg = p/{eB), i.e. to its momentum. 

For a chosen shock velocity and the magnetic field con- 
figuration we inject particles in the shock at some initial 
momentum po and follow their phase-space trajectories. We 
assume the constant particle injection to continue in time 
after the initial time to = 0. As some particles escape from 
the acceleration process by crossing the escape boundary 
placed far behind the shock we use the trajectory splitting 
procedure to keep the total amount of particles involved 
in simulations constant (cf. Kirk & Schneider 1987b; Os- 
trowski 1991). Here we put the boundary at the distance 
()K2,n/U2 + 4rj,2 . We checked by simulations that any fur- 
ther increase of this distance does not infiuence the results 
in any noticeable way. For every shock crossing, the particle 
weight factor multiplied by the inverse of the particle ve- 
locity normal to the shock (= particle density) is added to 
the respective time and momentum bin of the spectrum, as 
measured in the shock normal rest frame. As one considers 
a continuous injection in all instants after to, in order to 
obtain the particle spectrum at some time tj > to, one has 
to add to particle density in a bin pi at tj the densities in 
this momentum bin for all the earlier times. The resulting 
particle spectra are represented as power-law functions with 
the squared exponential cut-off in momentum 

f{p,t) = Ap-''e-(^y . (3.4) 

In this formula three parameters are to be fitted: the nor- 
malization constant A, the spectral index for the stationary 
solution a, and the momentum cut-off Pc (Fig. 2; for details 
of the fitting procedure see Appendix A). 

In the simulations, due to our, proportional to momen- 
tum, scaling of the respective quantities, the derived accel- 
eration time scale (2.3) must be also proportional to p, and 
thus to rg (p) /c. Therefore, this time scale measured in units 
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Figure 2. Particle spectra at a sequence of time instants. Indi- 
vidual simulation points represent the particle number (weight) 
AN per logarithmic momentum bin A log p. With full lines we 
present the respective fits (3.4). 



of rg{pc)/c (or re(pc)/c) is momentum independout and can 
be easily scaled to any momentum. The parameter gives 
the value of the acceleration time scale in units of re,i/c, 
T^cl ~ Tr re,i/c. The value of T^^' - at a particular time U 
is derived from the respective values of pc,i'- 

rp(c) _ Pc,i /„ 

acc,i pf. i—pc i — 1 ' V ■ / 

ti-ti-1 

where we consider the advanced phase of acceleration 
{pc,i S> po). As in our simulations pc oc t the condition 
{Pc,i — Pc,i-i) /Pci <C 1 is not required to hold in equation 
(3.5). Therefore, with all scales proportional to the parti- 
cle momentum, the formula (3.5) reduces to T^°^ ^ = U and 
the parameter Tr tends to a constant (Fig. 3). The exten- 
sion of the simulated spectra over several decades in particle 
energy allows to avoid problems with the initial conditions 
and decrease the relative error of the derived time scale by 
averaging over a larger number of instantaneous Tacc,i- 



4 THE ACCELERATION TIME SCALE IN 
RELATIVISTIC SHOCK WAVES 

For a given relativistic shock velocity particle anisotropy in 
the shock depends on the mean magnetic field inclination 
to the shock normal and the form of turbulent field. Below, 
we describe the results of simulations performed in order to 
understand the time dependence of the acceleration process 
in various conditions. In order to do that, we consider shock 
waves propagating with velocities Ui = 0.3, 0.5, 0.7 and 0.9 
of the velocity of light, and the magnetic field inclinations 
including the quasi-parallel, oblique sub-luminal and oblique 
super-luminal configurations. In all these cases we investi- 
gate the role of varying magnitude of turbulence character- 
ized here by the value of At or by the ratio of the diffusion co- 
efficient across the mean field and that along the field, k,± 
. The relation between these parameters for AQ = 10° is 



f c) 

Figure 3. An example of simulated values of T^cc in units of 
Te.i/c in a single run of the code. The simulation time t is given 
in units of rg^i{p = po)/c in the shock normal rest frame. We fit 
the final acceleration time scale only to the points at the advanced 
phase of acceleration (dashed line) . The dispersion of these points 
defines the fitting error. 



presented at Fig. 4, where - at At > 0.01 - the presented 
relation has the power-law form k±/k|| = 6.3 • 10"^ (At)^ 

4.1 Parallel shocks 

The most simple case for discussion of the first-order Fermi 
acceleration is a shock wave with parallel configuration of 
the mean magnetic field. As an example we consider the 
shock with negligible field inclination tpi = 1°. For such a 
shock, the present simulations do confirm the expected relar 
tion of decreasing the acceleration time scale with increasing 
the shock velocity and the amplitude of trajectory pertur- 
bations (Fig. 5) . One should note at the upper panel of the 
figure that for short At the presented time scales decrease 
more and more slowly with decreasing At . It is due to the 
fact that starting from some value of At we reach condi- 
tions of nearly isotropic diffusion, « k±, and further de- 
creasing of the time delay between scatterings decreases the 
acceleration time in much the same proportion as the time 
unit re,i/c used to measure it (cf. Equ. 3.3, Fig. 4). In the 
lower panel of Fig. 5 the diffusive time scale To is used as 
the time unit. The minute differences between the succes- 
sive curves reflect the statistical fluctuations arising during 
simulations. Without such fluctuations all curves should co- 
incide. The one sigma fit errors of T^cl are indicated near 
the respective points. One should note that for increasing 
the shock velocity the acceleration time scale decreases with 
respect to the diffusive time scale. 

4.2 Variation of T^cc with magnetic field 
inclination 

In order to compare the acceleration time scales for differ- 
ent magnetic field inclinations ipi we performed simulations 
assuming a constant scattering parameters upstream and 
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Figure 4. The relation between the scattering parameter At and 
the respective ratio of diffusion coefficients k±/kh. 
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Figure 5. The values of tIic\ in parallel shock waves {ipi = 1°) 
in units of a.) re,i/c and b.) Tq versus the shock velocity Ui. 
The values resulting from simulations are given for Ui = 0.3, 0.5, 
0.7 and 0.9 with the scattering amplitude parameter At near the 
respective results at the upper panel. 



downstream, yielding the same ratio of «;x/k|| in these re- 
gions. However, due to shock compression the particle gy- 
ration period is shorter downstream than upstream, in pro- 
portion to the mean magnetic field compression (3.1). At 
Fig. 6 wc present the values of the acceleration time scale 
derived in such conditions at different tpi. For super- luminal 
shocks the results are presented for the cases allowing for 
particle power-law energy spectra, i.e. when the cross-field 
diffusion is sufficiently effective. Actually, the spectra with 
inclinations a < 10.0 arc only included. Wc consider the 
following values of the magnetic field inclination: tpi = 1°, 
25.8°, 45.6°, 60°, 72.5°, 84.3° and 89°. The first one is for a 
parallel shock, the last two ones arc for perpendicular super- 
luminal shock with all velocities Ui. The intermediate values 
define luminal shocks ([/i/cosV'i = 1-0) at the successive 
velocities considered, respectively Ui = 0.9, 0.7, 0.5 and 0.3. 

In general, the acceleration time scale decreases with 
increasing field inclination, reaching in some cases the val- 
ues comparable, or even smaller than the particle upstream 
gyroperiod (6.28 in our units of re,i/c). The trend can be 
reversed for intermediate wave amplitudes when the mag- 
netic field configuration changes into the luminal and super- 
luminal one. Such changes arc accompanied with the steep- 
ening of the spectrum (see below). The acceleration rate 
at different scattering amplitudes changes with tpi in a way 
that at different inclinations the minimum acceleration times 
occur at different perturbation amplitudes (different At). 

An important feature of the acceleration process in rela- 
tivistic shocks should be mentioned at this point. The varia- 
tions of Ticc in oblique shocks are accompanied by changes of 
the particle spectrum inclination (cf. Kirk & Heavens 1989; 
Ostrowski 1991). At Fig. 7, the curves at {Ticl, a) plane 
represent the results for decreasing the scattering amplitude 
expressed with parameter At, and joined with lines for the 
same magnetic field inclination ipi . For parallel shocks the 
changes in Tacl do not lead to any variation of the spec- 
tral index. However, for oblique sub-luminal {ipi = 25.8°, 
45.6°) and luminal (^i = 60°) shocks a non-monotonic be- 
haviour is seen. The trend in changing Tacl and a observed 
at smaller perturbation amplitudes (larger At) is reversed at 
larger amplitudes, when the substantial cross-field diffusion 
is possible. 



4.3 Variation of Tacl with varying turbulence level 

In parallel shocks the acceleration time scale reduces with 
the increased turbulence level in its neighbourhood. This 

phenomenon, well known for non-relativistic shocks (cf. La- 
gage & Cesarsky 1983), is confirmed here for relativistic 
shock velocities (Fig. 5). In general, there are two main 
reasons for this change. The first one is a simple reduction 
of the diffusion time of particles outside the shock due to 
shorter intervals between scatterings, analogous to the de- 
crease observed in non-relativistic shocks. However, the in- 
creased amount of scattering influences also the acceleration 
process due to changing (decreasing) the particle anisotropy 
at the shock and, thus, modifying the mean energy gain 
of particles interacting with the shock discontinuity. Addi- 
tionally, in oblique shocks the upstream-downstream trans- 
mission probability may increase. One should note that the 
present approach is not able to describe fully the effect of de- 
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Figure 6. The value of T^^c in units of re.i/c (upper panels) and To (lower panels) versus the magnetic field inclination i/ii. The values 
resulting from simulations are presented for Ui = 0.3 and 0.9, with values of the parameter At given near the respective results. 



creasing anisotropy with the small amplitude random scat- 
tering model applied. It is due to the fact that correlations 
between the successive modifications of a trajectory (a se- 
quence of small angle scattering acts in this paper) in a 
single MHD wave cannot be accurately modelled within the 
simplified approach used. A more exact approach requires 
integration along the particle trajectories in realistic config- 
urations of the magnetic field. However, the comparison of 
the present simplified method to the one involving such an 
integration shows a reasonably good agreement (Ostrowski 
1993) suggesting that averaging over realistic trajectories is 
equivalent in some way to such averaging within our random 
scattering approach. 



side the shock is related to the normal diffusion coefficientHJ 
Hn {iin,i ~ _i cos^ + K.± J siTi'^ tpi , 1 = 1, 2), the increas- 
ing K± will lead - for large scattering amplitudes - to longer 
Tacc in units of [ve/c] . In the units of To the acceleration 
time depends only weakly on the turbulence level and shows 
a small maximum for the minimum at the presented figure. 
For super- luminal shocks one can note absence of data points 
corresponding to low turbulence levels, where the power-law 
spectrum cannot be formed or it is extremely steep. In these 
excluded cases, the upstream population of energetic parti- 
cles is only compressed at the shock with the characteristic 
upstream time of ~ r^^i/Ui (cf. Begelman & Kirk 1990; Os- 
trowski 1993). 



In shocks with oblique magnetic fields a non-monotonic 
change of the acceleration time scale with the amount of 
scattering along the particle trajectory is observed (Fig. 8, 
see also Fig. 6; cf. Ostrowski 1991, 1994b for the spectral in- 
dex). Increasing the amount of turbulence up to some critical 
amplitude decreases the diffusion time along the magnetic 
field and thus Tacl- However, as the mean diffusion time out- 



t One should note that for the relativistic shocks, due to particle 
anisotropy, the respective relation may be not so simple as that 
given in equation (2.2) for non-relativistic shocks. 
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Figure 7. The values of T^cc in units of r^^^/c at different inclina- 
tions versus the particle spectral index a. The values resulting 
from simulations are given for Ui = 0.5 for five values of the angle 
tpi given near the respective results. The maximum value of At 
is given at the end of each curve and it monotonously decreases 
along the curve. 
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Figure 8. The acceleration time T^cc versus the level of particle 
scattering for shocks with velocity Ui = 0.3 (dashed lines) 

and 0.5 (full lines). We present results for three values of the 
magnetic field inclination: a.) parallel shock (i/ii = 1°), b.) a sub- 
luminal shock with i/ii = 45.6° and c.) a super-luminal shock with 
ipi = 89°. 



4.4 Large angle scattering (LAS) versus pitch 
angle diffusion (PAD) 

Some earlier derivations of the acceleration time scale were 
based on the numerical simulations involving particle scat- 
tering at point like scattering centres isotropizing the par- 
ticle momentum at each scattering. This approach does not 
provide a proper description for the acceleration processes 
in shock waves moving with velocities comparable to the 
particle velocity because it removes particle anisotropy and 
changes the factors related to it. Moreover, against argu- 



ments presented in some papers, such scattering pattern 
can not be realized in turbulent magnetic fields near rel- 
ativistic shocks, where most particles active in the accel- 
eration process are able do diffuse only a short distance, 
below a few particle gyroradii off the shockla. Such distances 
are often insufficient to allow for big particle pitch-angle 
changes occuring with the point-like scattering centres which 
isotropize particle momentum at each scattering. In shocks 
with oblique magnetic fields such scattering pattern can sub- 
stantially change the shape of the accelerated particle spec- 
trum with respect to the PAD model. Additionally, as an 
individual particle interaction with the shock can involve 
a few revolutions along the magnetic field, the usually as- 
sumed adiabatic invariant conservation, p\/B — const, can- 
not be valid for short inter-scattering intervals. 

Before proceeding with the results let us remark a fur- 
ther problem arising within any model involving discrete 
particle scattering acts. As discussed in Section 3 the pres- 
ence of scattering is equivalent to the presence of magnetic 
field perturbations. As a result the effective field is larger 
than the uniform background field. The presented estimate 
of the this field (Equ. 3.3) is valid for small angle momen- 
tum scattering. The amount of energy in magnetic turbu- 
lence with the waves shorter than cAt is required to be 
small because the presented estimate assumes the particle 
momentum perturbation in At occurrs on the uniform ef- 
fective perturbing field. To compare the scattering processes 
with different At one have to neglect the unknown factor of 
the ratio of the averaged actual magnetic field to the esti- 
mated value (like the one in Equ. 3.3). Let us note that this 
factor, as well as the notion of the effective field were not 
considered in the earlier papers. In derivation of the results 
presented below we considered either the Equ. (3.3) to be 
formally valid for both, PAD and LAS scattering models, or 
we presented data with the use of units defined by the parti- 
cle gyroradius in the uniform magnetic field component Bq. 
The later case is provided to be compared with the earlier 
results by other authors. 

In order to have the comparison of the LAS versus PAD 
models meaningful we simulate the acceleration process for 
LAS with the described in Section 3 procedure, but applying 
the isotropic scattering of particle momentum. Instead of 
our An = 10° for the pitch angle diffusion we use Af2 = 
180° for the large angle scattering. The times At between 
successive scatterings are chosen in a way to have the same 
spatial diffusion coefficients along the magnetic field for both 
models, as measured in the units related to the effective 
field component (re and re/c). With the formulae derived 
in simulations (^f'*'" = 41.7Atio° for PAD and Ky"*'^ = 
0.174Ati8o° for LAS) given in the units based on Bo {rg 
and Tg/c), and with the equation (3.3) for the effective field 
Be one obtaines the required relation of the scattering times: 
AtisQo = [57611Atioo + 775.4]^/^. Let us note that as our 
procedure derives the actual particle trajectories (not only 
the gyration centre motion) any random scattering process 
results in the cross-field diffusion and breaking the strict 
p\/ B conservation. 

t However, for the non-relativistic shock velocity and particles 
much above the injection energy such approximation can be safely 
used (cf. Jones & Ellison 1991). 
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Figure 9. Comparison of the acceleration times derived with our 
pitch angle diffusion procedure, T^J?^ = T^'c'ci^^ = 10°), and 
the large angle scattering model, T^fO = rici(An = 180°). At 
the figure the points for the shock velocity Ui = 0.5 and two 
magnetic field inclinations, *i = 0° and 45.6°, are joined with 
lines for successive values of At . Pull (dashed) lines are for the 
data in the units defined by Be (-Bo), respectively. The points for 
maximum At (= 6.28 in PAD) are indicated near the respective 
curves. 
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Figure 10. Comparison of the spectral indices derived with our 
pitch angle diffusion procedure, a^" s a(Af2 = 10°), and the 
large angle scattering model, a^*" = a{ACl = 180°). The data 
are presented for the shock wave with Ui = 0.5 and the magnetic 
field inclinations 0° and 45.6°. At the figure the data are plotted 
in the same way as at Fig. 9. 



The discussion in this sub-section is to clarify the re- 
lation of our present simulations to some earlier ones pro- 
viding sufficient details on the applied numerical approach. 
We performed simulations with the use of the LAS model 
for Ui — 0.3, 0.5, 0.7 and 0.9, and two inclination angles, 
ipi = 0.1° and 45.6°. The results for Ui = 0.5 are presented 
at Fig-s 9 and 10. The comparison to the earlier results done 
for the largest At values, where the effective field approx- 



Figure 11. The ratio of the PAD simulated time scale Tac^ 
to the 'theoretical' non-relativistic time scale T^cc°^ given in 
Equ. 2.2 for different inclinations of the magnetic field 'I'l. The 
results for Ui = 0.3 and 0.5 are presented for the case of negligible 
cross-field diffusion (At = 6.28 in the PAD model). 



imately equals the background field, shows a general qual- 
itative agreement with the results of Ellison et al. (1990) 
and Naito & Takahara (1995). For parallel shocks the accel- 
eration time scale for LAS is only slightly longer than the 
PAD one for all considered shock velocities. The ratio of this 
time to the time scale estimate based on the equation (2.1) 
for non-rclativistic shocks yields values decreasing from ap- 
proximately 0.6 at Ui = 0.3 up to 0.17 at Ui — 0.9. At high 
velocities the ratio is much smaller (a factor of 2 - 3) than 
the ones derived by Ellison et al. (1990). However, there ex- 
ist a substantial difference of numerical procedures applied 
here and in the mentioned paper, and we are not able to 
point out any single reason for the noted discrepancy. As 
mentioned previously our time scale defined with the use of 
the cut-off energy is shorter than the steady-state estimate 
derived by Ellison ct al. Additionally, in the mentioned paper 
the use of analytic solutions of the spatial diffusion equation 
for the case of strongly relativistic flow can be questioned. 
For oblique shocks we observe analogous reduction of the 
acceleration time scale as that reported by Naito & Taka- 
hara (1995), with the PAD model allowing for a little more 
rapid acceleration than the LAS model. Of course this agree- 
ment is broken for short At, where the cross field diffusion 
can not be neglected and the particle magnetic momentum 
is not conserved at interactions with the shock. At Fig. 10 
one can compare the spectral indices derived with the use of 
both approaches. As a final remark we would like to point 
out that the only independent check for any numerical pro- 
cedure is provided by comparying of the derived spectral 
indices in the limiting cases of SB « B and n± ~ Ky to 
the analytic results (e.g. Heavens & Drury (1988) for parallel 
shocks and Kirk & Heavens (1989) for oblique shocks). This 
comparison for our PAD method is presented in Ostrowski 
(1991). 
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5 SUMMARY AND FINAL REMARKS 

In the present paper we consider the acceleration time scale 
in relativistic shock waves. We argue that the numerical ap- 
proach based on the model involving particle isotropic scat- 
tering neglects or underestimates a very significant factor 
controlling the acceleration process - the particle anisotropy. 
We also note that such a scattering pattern cannot be physi- 
cally realized when particle trajectory perturbations are due 
to MHD waves. It is why the derivations of the present 
work are based on the small-angle momentum scattering ap- 
proach. We note that the code allows for a satisfying descrip- 
tion of particle transport in the presence of large amplitude 
magnetic field perturbations as well. We demonstrate the ex- 
istence of correlation between particle energy gains and its 
diffusive times. The analogous correlation is expected for the 
probability of particle escape downstream the shock. There- 
fore, for defining the acceleration time scale we use the rate 
of change of the spectrum cut-off momentum which acom- 
modate all such correlations. We performed Monte Carlo 
simulations for shock waves with parallel and oblique (ei- 
ther, sub-luminal and super-luminal) magnetic field config- 
urations with different amounts of scattering along parti- 
cle trajectories. Field perturbations with amplitudes ranging 

from very small ones up to SB ~ B arc considered. In par- 
te) 

allel shocks Tacc diminishes with the growing perturbation 
amplitude and the shock velocity. However, it is approxi- 
mately constant for increasing turbulence level if we use the 
respective diffusive time scale as the time unit. Another fea- 
ture discovered in oblique shocks is that due to the cross-field 
diffusion T^cl can change with SB in a non-monotonic way. 
The acceleration process leading to the power-law particle 
spectrum in a super-luminal shock is possible only in the 
presence of large amplitude turbulence. Then, the shorter 
acceleration times occur when the perturbations' amplitudes 
are smaller and the respective spectra steeper. We discussed 
the coupling between the acceleration time scale and the par- 
ticle spectral index in oblique shock waves with various field 
inclinations and revealed a possibility for non-monotonic re- 
lations of these quantities. 

The shortest acceleration time scales seen in the Simula^ 
tions are below the particle gyroperiod upstream the shock. 
These times do not require the ultra-relativistic shock veloci- 
ties, but may occur in mildly relativistic ones with the quasi- 
perpendicular magnetic field configuration. One should note 
that - due to the larger magnetic field downstream the shock 
- in this short time the particle trajectory can follow a few 
revolutions near the shock with only a short section of each 
one penetrating the upstream region. 

The presented estimates of the acceleration time scale 
provide an interesting possibility for modelling shock waves 
in the conditions where the electron spectrum cut-off en- 
ergy is determined by the balance of gains and losses. If 
one is able to derive the respective acceleration rate from 
the knowledge of the energy loss process, and the particle 
spectral index is also known, then both these values provide 
constraints for the acceleration process which could be fur- 
ther used to reduce the parameter space available for the 
considered shock wave (cf Fig. 7). 

The work of MO was partly done in Max-Planck- 
Institute fiir Radioastronomie in Bonn. He is grateful 
to Prof. R. Wielebinski for a kind invitation. The pre- 



sented computations were done on the CONVEX Exemplar 
SPPlOOO/XA-16 in ACK 'CYFRONET' in Krakow. The 
work was supported by the grant PB 1117/P3/94/06 from 
the Komitet Badan Naukowych. 



APPENDIX A. FITTING THE PARAMETERS 

A, a AND Pc 

Any simulated spectrum was evolving in time by increasing 
the width of its power-law section and, thus, the best fit of 
this power-law was possible with the use of the final spec- 
trum at maximum time. Therefore, in the simulations we 
used the last spectrum to fit parameters of the power-law 
normalization and the spectral index, A and a. Next, for 
any earlier spectrum, these parameters were assumed to be 
constant and we were fitting only the cut-off momentum, 
Pc- For each fit we used 20 last points of the spectrum pre- 
ceding the point where particle density fell below 0.16 of 
Ap~". The number of 0.16 was chosen experimentally as to 
obtain the best fits to the cut-off region of the spectrum. As 
the distribution (3.4) represents only an approximation to 
the actual particle distribution, there was no reason to use 
points corresponding to lower densities, of lesser statistical 
significance. 
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